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ABSTRACT 

Algol is a triple stellar system consisting of a close semi-detached binary orbited by a third 
object. Due to the disputed spatial orientation of the close pair, the third body perturbation 
of this pair is the subject of much research. In this study, we determine the spatial orientation 
of the close pair orbital plane using the CHARA Array, a six-element optical/IR interferometer 
located on Mount Wilson, and state-of-the-art e-EVN interferometric techniques. We find the 
longitude of the line of nodes for the close pair is fii — 48° ± 2° and the mutual inclination of 
the orbital planes of the close and the wide pairs is 95° ± 3°. This latter value differs by 5° from 
the formerly known 100° which would imply a very fast inclination variation of the system, not 
borne out by the photometric observations. We also investigated the dynamics of the system 
with numerical integration of the equations of motions using our result as an initial condition. 
We found large variations in the inclination of the close pair (its amplitude ~ 170°) with a period 
of about 20 millenia. This result is in good agreement with the photometrically observed change 
of amplitude in Algol's primary minimum. 



1. Introduction 

There are about 1000 triple stellar systems 
known in the Galaxy, many of which consist of 
a close eclipsing pair and a distant third ob- 
ject orbiting around the close pair (Batten 1973, 
Tokovinin 1997). Algol is probably the most well- 
known of such systems. 

Algol consists of a semi-detached eclipsing bi- 
nary with an orbital period of 2.87 days (B8V 
+ K2IV) with an FlIV spectral type star revolv- 
ing around the binary every 680 days (discovered 
by radial velocity measurements, Curtiss 1908). 
Early interferometric observations were unable to 
resolve the system (Merrill 1922), but the third 
component was succesfully observed by speckle in- 
terferometry (Gezari, Labeyrie and Stachnik 1972, 
Blazit et al. 1977, McAlister 1977, 1979) and its 
orbit was precisely determined by Bonneau (1979). 
This result was refined by using the Mark III op- 
tical stellar interferometer (Pan et al. 1993). 

In the radio regime, Lestrade et al. (1993) de- 
tected positional displacement during the orbital 
revolution of the AB pair using the VLBI tech- 
nique, and identified the K-subgiant as the source 
of radio emission. The orbital elements of the close 
and the wide pairs determined from all these ob- 
servations are listed in Table [T] 

The light minima - mainly primary - of Algol 
were extensively observed in the last two centuries. 
There was only a very small change in the eclipse 
depth during this time. This led Soderhjelm (1975, 
1980) to the theoretical conclusion that the mutual 
inclination of the orbital planes of the close and 
the wide pair systems should not be larger than 
11° and likely they are coplanar because both the 
shape and depth of the light minima should have 
noticeably changed otherwise. This theoreti- 
cal result was seemingly in good correspon- 
dence with the inclination data deduced 
both for the close and the wide orbits, i.e. 
ii = 82?3 ± 0?2, and ^2 = 83° ± 2° (Soderhjelm 
1980) (see also Figure 1), respectively, from 
which Soderhjelm (1980) stated the exact 
coplanaritj0. Based on the earliest speckle 
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At this point we should take a clear distinction be- 
tween the different kind of orbital elements which 
are mentioned in this paper, and describe our nota- 



interferometric measurements, Soderhjelm 
(1980) calculated ^2 = 132° ± 2° for the 
wide orbit, and therefore he expected f^i = 
132° ± 4° for the node of the eclipsing pair. 
Later Pan et al. (1993) determined the 
astrometric orbit of the third component 
and found that both the longitude of the 
ascending node {fl2) and the argument of 
the periastron {uj2) of the wide orbit practi- 
cally differ by 180° from the previously ac- 
cepted values. In the case of an isolated 
two-body astrometric orbit, 180° discrep- 
ancy is non-problematic, because geomet- 
rically it means the reflection of the orbital 
plane onto the plane of the sky for which 
transformation the astrometric coordinates 
are invariants. Nevertheless, in a triple 
body system this results in different spatial 
configuration of the orbital planes, i.e. it 
modifies the mutual inclination fundamen- 
tally. However, the polarimetric measure- 
ments of Rudy (1979) yielded a contradic- 
tory result, suggesting that fii — 47° ± 7° 
which would imply a perpendicular rather 
than coplanar configuration. Nevertheless, 
Soderhjelm (1980) suggested, that perhaps 
the nature of the polarization mechanism 
was not understood correctly. Note that 
Rudy (1979) resolved also the inclination 
ambiguity, i.e. he determined that the 
angular momentum of the binary directed 



tion system. The optical interferometric measure- 
ments give information about the relative motion 
of one component to the other. From the CHARA 
measurements we get information about the rela- 
tive orbit of Algol B around Algol A. The orbital 
elements refer to this relative orbit denoted by sub- 
script 1. Similarly, the earlier astrometric measure- 
ments of the third star give its relative orbit to the 
close binary. (More strictly speaking, to its photo- 
center.) These elements are denoted by subscript 
2. VLBI measurements give the motion of Algol 
B component in the sky, i.e. after the use of the 
necessary corrections we get the orbital elements 
of the secondary's orbit around the center of mass 
of the binary. These elements are denoted by sub- 
script B- Finally, in order to carry out some of 
the aforementioned corrections for calculating the 
orbital motion of Algol B, we need the orbital ele- 
ments of the close binary in its revolution around 
the center of mass of the whole triple system. These 
orbital elements are denoted by Nevertheless, 
the elements of these latter two orbits will be only 
used when necessary. 
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away from the observer, and, consequently, 
the system incUnation (u) should be less 
than 90°. 

This nearly perpendicular configuration was 
supported by other measurements: Lestrade et 
al. (1993) found Qi = 52° ± 5° for the close 
pair in good agreement with the polarimetric mea- 
surements of Rudy (1979). If this value of Oi is 
correct, then the mutual inclination is about 100° 
(Kiseleva et al. 1998), and the two orbital planes 
are nearly but not exactly perpendicular to each 
other. This value for the mutual inclination of the 
system has been widely accepted since then. How- 
ever, we propose that this mutual inclination 
value cannot be correct due to dynamical consid- 
erations which is described hereafter. 

It is well-known that in a hierarchical triple 
stellar system, the orbital planes of the close and 
wide pairs are subject to precessional motion, in 
such a way that the normals of the orbital planes 
move on a conical surface around the normal of 
the invariable plane. (Here we omit the effect of 
stellar rotation which is insignificant in an ordi- 
nary triple system.) In the case of a hierachical 
triple system where the invariable plane almost 
coincides with the wider orbital plane, the preces- 
sion cone angle is close to the mutual inclination. 
Consequently, in the case of the present system, we 
would get an almost 160°-amplitude variation in 
the observable inclination during the approximate 
period given by Eq. (27) of Soderhjelm (1975) (as- 
suming the present approximation is valid, when 
the mutual inclination tends to 90° , the precession 
period tends to infinity). Here we refer to Fig. 4 of 
Borkovits et al. (2004) which clearly shows that in 
the case of the aforementioned configuration the 
observable inclination of the close binary would 
have changed by approximately 3° in the last cen- 
tury which evidently contradicts the observations. 
In this case the eclipses would disappear within a 
few centuries. This was already observed in some 
eclipsing binaries, like in SS Lacertae or V907 Sco 
(see e.g. Eggleton & Kiseleva-Eggleton 2001) but 
not in Algol. 

In summary, the polarimetric and the inter- 
ferometric observations contradict the coplanar- 
configuration, but - since the mutual inclination is 
far from exact perpendicularity - the latter is not 
in agreement with the observed tiny change in the 
minima depth. A closer approximate perpendicu- 



larity of the two orbital planes would mean that 
the period of orbital precession becomes so large 
that the inclination variation (and consequently 
the depth variation of the minimum) of the close 
pair remains unobservable for a long time, consis- 
tent with the observations. 

The aim of this study was to constrain the mu- 
tual inclination of the system better, requiring the 
measurement of the longitude of the node for the 
close pair. The other orbital elements are well- 
known from spectroscopic or photometric data, 
but there is a controversy in the value of fii. Be- 
cause the expected apparent size of the close bi- 
nary semi-major axis is of the order of 2 milliarc- 
seconds (mas), we carried out optical and radio 
interferometry measurements. As we will show, 
optical and radio interferometry are complemen- 
tary techniques. Combining these two, we will 
show that it is possible to resolve the ambiguity 
in the geometry of the system, and better assess 
the accuracy of our measurements. 

2. Observations and data reduction 

2.1. CHARA Observations 

The CHARA Array is an optical/ncar-IR in- 
terferometer array consisting of six Im telescopes. 
The array is described in detail in ten Brumme- 
laar et al. (2005). A detailed overview and further 
references about the observables and the theory 
of optical interferometry can be found in Haniff 
(2007). 

We observed Algol on three nights (2, 3 and 4 
December, 2006) in the Ks band (the effective 
wavelength was 2.133 /im). Much of the sec- 
ond night was lost due to high winds and dusty 
conditions. 

Iota Persei and Theta Persci served as calibra- 
tion stars. The observations of the target and the 
two calibrators were organized into a sequence and 
the measurenicnits on calibrators generally brack- 
eted the target observations. We have 12 data 
points of Iota Persei, 12 data points of Theta Per- 
sei and 23 data points of Algol itself 

Each data point was calculated from a number 
of scans. Each scan measured the intensity vari- 
ations as a function of the path delay. About 300 
scans were collected within five minutes for one 
data point, the first 22 of these were obtained on 
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the targets (Algol or one of the calibrator stars), 
then 21 scans were done for measuring the back- 
ground while the shutter was closed, then more 
than 200 other scans were obtained on the targets 
again and finally 67 further scans for measuring 
the background again. 

To reduce the data, we used the recipe of McAl- 
ister (2002). This consisted of the following steps: 
First a low-pass filter was applied to remove the 
atmospheric noise. Then the bias was subtracted 
and the scans were normalized to unity. As a 
next step, the scans measured by the two chan- 
nels were subtracted from each other (for details, 
see ten Brummelaar et al. 2005 and McAlistcr 
2002). This was further processed by applying 
a high-frequency filter to reduce noise. Discrete 
Fourier-transforms of the scans were calculated 
and a template was computed. For this new tem- 
plate, we used the full amplitude for the frequen- 
cies ±25/cycle around maximum frequency and 
20% of the amplitude for the other frequencies 
(for more details see McAlister 2002). From 
this we could calculate the maximum deviation 
of the template from zero which yielded an esti- 
mation of the visibility value. Removing outlier 
values, we averaged the remaining values which 
yielded the uncalibrated visibility of a particular 
point. The error was estimated as the standard 
deviations of the visibility values of the more than 
200 scans of the point. 

The measured visibilities of the calibrators were 
linearly interpolated for the times of the Algol ob- 
servations. A comparison of the true and mea- 
sured visibilities of the calibrators yielded a factor 
which converted the measured target visibilities to 
true ones. The true visibilities of the two calibra- 
tors were estimated as follows: 

The uniform disc (UD) angular diameter of Iota 
Persei is 1.21 ± 0.06 mas according to the 'Cat- 
alog of High Angular Resolution Measurements' 
(Richichi et al. 2005). The true diameter of 
the other calibrator star Theta Persei is 1.201 ± 
O.O15i?0 which was determined by the fit of its 
spectrum (Valenti and Fisher 2005). Since its 
HIPPARCOS parallax is known one can easily cal- 
culate its true UD angular diameter to be 0.995 ± 
0.01 mas. 

The limb-darkened (LD) angular diameter is 
larger than the UD diameter. There exists a sim- 
ple relationship between them (Hanbury Brown et 



al. 1974): 

Old _ / 1 - ux/3 

dvD^^ I- 7ua/15 ^ ' 

The limb-darkening coefficients for both compo- 
nents were taken from the tables of van Hamme 
(1993). These coefficients are a function of sur- 
face gravity and effective temperature which them- 
selves were estimated from the known spectral 
type of the cahbrators. Then Eq. ([1]) yielded the 
corrections which increase the UD angular diame- 
ters by a few percent only. These corrected values 
were used to calibrate the visibilities. 

The telescope combinations, epoch of observa- 
tions, baselines, uncalibrated visibilities and their 
errors can be found in Table [21 The UV-coverage 
can be seen in Fig. [5] 

2.2. e-VLBI Observations 

We observed Algol with a subset of the Euro- 
pean VLBI Network (EVN) betv^^een 16:37 and 
1:19 UT on 14-15 December 2006 at 5 GHz. 
These observations were carried out using the 
e-VLBI technique, where the telescopes stream 
the data to the central data processor (JIVE, 
Dwingeloo, the Netherlands) instead of recording. 
Using e-VLBI for the observations was not fun- 
damental for our measurements, but we took the 
opportunity of an advertised e-VLBI run close 
in time to the CHARA observations, outside the 
normal EVN observing session. The participat- 
ing telescopes were Cambridge and Jodrell Bank 
(UK), Medicina (Italy), Onsala (Sweden), Torun 
(Poland) and the Westerbork phased array (the 
Netherlands). The data rate per telescope was 
256 Mbps, which resulted in 4x8 MHz subbands 
in both LCP and RCP polarizations using 2-bit 
sampling. The correlation averaging time was 2 
seconds, and we used 32 delay steps (lags). Initial 
clock searching was carried out before the experi- 
ment using the fringe-finder source 3C345. Algol 
was phase-referenced (Beasley & Conway 1995) 
to 0309-1-411 in 3-5-3 minute cycles. Additional 
scans were scheduled on 3C84 for real-time fringe 
monitoring, and for D-term calibration. We used 
3C138 to calibrate the Westerbork synthesis array 
amplitudes and polarization. 

Post-processing was done using the US Na- 
tional Radio Astronomy Observatory (NRAO) 
AIPS package (Diamond 1995). The amplitudes 
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were calibrated using the known antenna gain- 
curves and the measured system temperatures. 
The data were fringe-fitted, bandpass caUbrated, 
and then polarization calibrated. We corrected 
for the polarization leakage D-terms and fringe- 
fitted the cross-hand data, after which the data 
were averaged in frequency in each subband. Be- 
sides the standard procedure, we used the WSRT 
synthesis array measurements on 0309-1-411 to 
obtain a more accurate VLBI flux scale. The 
phase-reference source showed a low level of cir- 
cular polarization (fractional CP ~ 0.28%). The 
left and right-handed VLBI gains were separately 
adjusted in accordance with the WSRT mea- 
surement. Imaging was carried out in Difmap 
(Shepherd et al. 1994). The snap-shot images 
(from about 45 minutes data each) were made 
by Fourier-transforming the observed visibilities, 
no self-calibration was applied. We fit circular 
Gaussian model components to the ww-data 
in Difmap. Initially one component v^ras fit 
in each snapshot. Then, the size and flux 
of the component v^ras fixed, and we let the 
position vary for each 5-minutes Algol scan. 

The log of the observations as well as the cal- 
culated individual relative positions can be found 
in Tabled 

3. Data analysis and results 
3.1. Analysis of CHARA data 

According to the van Cittert-Zernike theorem, 
the amplitude of the visibility is the normalized 
Fourier-transform of the intensity distribution (for 
a comprehensive explanation, see Haniff 2007): 



V{u,v,t) = 



J J I{x,y,t) cos (2n^ 



AxAy 



J J I{x,y,t)dxdy 



(2) 

In this equation V is the true visibility at the (u, 
v) spatial frequencies, (x, y) are the corresponding 
sky coordinates, t is the time, / is the intensity at 
the {x,y) sky point and finally A = 2.133 /im is 
the effective wavelength of our observations (cor- 
responding to the Kg band). Note that the inten- 
sity distribution on the sky normally changes very 
slowly with time, but we had to include the time 
dependence in Eq. because of the rapid orbital 
revolution of the AB pair. Because of the power 
of CHARA, it is not enough to model the system 



as the combination of two point-like sources but 
we need to combine the pictures of two extended 
sources. 

To determine this intensity distribution we de- 
veloped a model which was very close to the one 
of Wilson and Devinney (1971) which is based on 
the Roche-model (Kopal 1978) but it was imple- 
mented in IDL to restore the surface intensities 
into a matrix and to calculate the sky-projected 
picture of the system. 

Since we have only about two dozen visibility 
measurements, we wanted to limit the number of 
free parameters, choosing just three: the angular 
size of the semi- major axis (in mas), the Ks surface 
brightness ratio, and the angle fii. All other pa- 
rameters were fixed according to the values given 
in Wilson et al. (1972) or in Kim (1989). 

The model outlined above was used to fit the 
data and a grid search was carried out on the free 
parameters: 

a) The surface brigthness ratio was stepped 
from to 1.4 with a step size of 0.05; 

b) rii was stepped from 0° to 180° with a step- 
size of 5° (note that the visibility amplitude 
does not change if we rotate the image by 
180°, which leaves a 180° ambiguity in the 
ascending node; this ambiguity in fii can 
be resolved with VLBI); and 

c) the angular size of the semi-major axis was 
stepped from 1.80 mas to 3.60 mas with a 
step-size of 0.005 mas (note that the ex- 
pected size was 2.50 mas). 

Nearly 353 000 models were calculated on this 
grid, and the minimum was found. Around the 
minimum, a new search was carried out with a 
finer grid and about 35 000 new models were com- 
puted again. Around the minimum, a polynomial 
fit yielded the final values and errors. The results 
are shown in Figs. [3] and ID 

The best solution we found is reported in Ta- 
ble SI 

3.2. Analysis of e-VLBI data 

The EVN measurements were carried out on 
14/15 December 2006 for almost 9 hours. As 
the orbital period of the eclipsing binary is P = 
2.8673 days, the orbital arc covered was near 13%. 
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During this period a secondary minimum occurred 
which was simultaneously observed photometri- 
cally by the 50-cm telescope at Piszkesteto Station 
of the Konkoly Observatory, Hungary (see Bfro et 
al., 2007). One of our goals was to observe possi- 
ble partial occultation of the radio source by the 
primary component and measure the change in 
the circular polarization properties of the source 
accordingly. Because Algol flared during the ob- 
servations (see Fig. [5]), this goal could not be 
fulfilled and will not be discussed further here. 

In contrast to the measurements of Lestrade et 
al. (1993) who observed the Algol on four differ- 
ent nights in near quadrature phases (i.e. around 
(j> = 0.25 and (f) = 0.75), our observations covered 
a small part of the orbit around (p ^ 0.5, when 
the arc projected onto the plane of the sky is the 
largest. (Naturally the case is the same at ~ 0). 
The advantage of this approach is that the mo- 
tion of the target can be detected in a few hours, 
and the measured positions are only slightly af- 
fected by the orbital motion of the AB-C pair, un- 
like the case when the data arc taken at different 
epochs. Because of the short observing time inter- 
val, the orbit of the AB pair itself is not well con- 
strained. Nevertheless, using a priori known 
values for most of the other orbital param- 
eters, one expects to find a relatively accu- 
rate value for the longitude of the ascending 
node. 

There are two limitations that must be men- 
tioned here, 1) short-term tropospheric and iono- 
spheric phase fluctuations which cannot be mod- 
elled well and limit the astrometric accuracy in 
short VLBI-measurements, and 2) the variable 
structure of the radio source (Mutel et al. 1998) - 
the radio emission is not coming from the surface 
of the K-subgiant, but likely from its active polar 
coronal region. Although Algol was unresolved 
with our array configuration, the source centroid 
position could have changed appreciably because 
of the bright flare during the run. These are the 
factors that must be taken into account in the in- 
terpretation of the final result. 

We calculated astrometric orbits as well as a 
simple linear fit on two different sets of observing 
data. First, hourly normal points were formed. 
Then we also calculated the orbit using 5 min- 
utes averages. The latter showed that some points 
with extremely large scatter can lead to false re- 



sult in the hourly normal points. Removing such 
outliers, we calculated our final solution from the 
5-minute-average data. Despite the shortness of 
our observing session (less than 9 hours) the posi- 
tional data were corrected for the annual parallax, 
proper motion and the revolution around the cen- 
tre of mass of the triple system. These corrections 
resulted in approximately 1° difference in the node 
position. For this latter correction we recalculated 
the third body orbit by the same code which was 
used for the binary orbit determination from our 
e-VLBI data. We used both the data sets of Bon- 
neau (1979) and Pan et al. (1993). From these 
(very similar results) we applied the orbital ele- 
ments obtained from Pan's data (see Table [5]) for 
the wide orbit correction. For the astrometric 
calculations, we used our own differential correc- 
tion code based on a Levenberg-Marquardt algo- 
rithm, tested against the data and results of Eich- 
horn and Xu (1990), and found it to be in excellent 
agreement. In our code the maximum num- 
ber of adjustable free parameters is nine: 
the position of the centre of mass {Xq, Yq), the 
orbital period (P), and the six usual orbital ele- 
ments (a, e, i, u!, il, and Mq). In the case of the 
close binary orbit determination, because of 
the circular orbit, instead of the argument of peri- 
astron (w), and the mean anomaly (M), their sum, 
i.e., the true longitude (the distance from the as- 
cending node in case of circular orbits) should be 
used. This was done in such a way, that uji was 
formally considered as zero, while (Afo)i = 90° 
was set for the mid- (secondary) eclipse moment, 
to = 2 454 084.360. The period and the inchnation 
were acquired from Kim (1989), while the semi- 
major axis of the secondary's orbit around the 
centre of mass of the binary {as) was calculated 
from Kim's data. First we adjusted three pa- 
rameters {Xq, Yq, fti), leaving the other six 
parameters fixed. Finally, the semi-major 
axis (as) was also adjusted, as a fourth pa- 
rameter. 

Using Gnuplot, we also performed a simple lin- 
ear fit to the data from the knowledge that in the 
vicinity of (j) ~ 0.5 the motion can be approxi- 
mated well with a line whose slope gives tanfii. 

Our solutions can be seen in Table O as well as 
in Fig. [6l 
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4. Discussion 

4.1. Discussion of the results 

First, we concentrate on the CHARA results. 
The true size of the semi-major axis of Algol is 
ai = 14.1 Rq (Kim 1989) while its angular size 
was measured by us to be ai = 2.28 ± 0.02 mas 
(see Table 13]). By dividing these two numbers one 
can find that the distance to Algol is 28.6±0.3 par- 
sec. The HIPPARCOS parallax yielded 28.46tg;^'^ 
parsec. The agreement is excellent. 

The determined surface brightness ratio (0.33) 
can be converted to a luminosity ratio by multi- 
plying with the ratio of surface area of the two 
stars, which is known from the light curve solu- 
tion (Wilson et al. 1972; Kim 1989). The resulting 
luminosity ratio is 0.43. This value is very close 
to the photometrically estimated 0.44 (Murad and 
Budding 1984). 

According to the CHARA results, the longitude 
of the node is Sli = 48° ± 2° with an ambiguity of 
180 degrees. Because the determined distance and 
luminosity ratios agree very well with earlier mea- 
surements obtained with other methods, we have 
confidence in our results. With the VLBl mea- 
surements (see below) we resolve the ±180° am- 
biguity and conclude that Qi = 48° ± 2°. This is 
in excellent agreement with the value determined 
from polarimetric measurements (fii = 47° ± 7°, 
Rudy 1979), indicating that polarimetry is an ef- 
ficient tool to determine the spatial orientation of 
the orbits. 

At this point we can determine the mutual in- 
clination i„i with the following formula: 

cosim = cosil C0SZ2 + sin?i sin 12 cos(51i — Q2) 

(3) 

The result is im = 95° ± 3° (the uncertainty re- 
fiects the uncertainties not only in but also in 
the other angular elements), confirming the con- 
clusion of Lestrade et al. (1993) that the two 
orbital planes are nearly perpendicular to each 
other. This value is, however, closer to the exact 
perpendicularity than the 100° given in Kiseleva et 
al. (1998) which was based on the measurements 
of Lestrade et al. (1993). Nevertheless, the exact 
perpendicularity is within the three sigma range. 

Regarding the e-VLBI measurements (see Ta- 
ble O, we obtained fii = 53° ± 826° and 
ni = 53° ± 239° for the node from the three 



and four adjusted parameter astrometric 
fits respectively, and fii = 52° ± 3° from the 
simple linear LSQ fitting. These values are close 
to those obtained previously. However, we note 
that in the case of the astrometric orbit fittings, 
the formal errors are extremely large. This natu- 
rally reflects the fact that our measurements cover 
only a very short fraction of the orbit, and espe- 
cially in that phase, where the expected astromet- 
ric orbit is almost a straight line. Consequently, 
without any a priori information, the orbit would 
be completely undeterminable. However, in this 
particular case, the longitude of the node itself is 
very well determined during this phase, as this 
is nothing other than the slope of the obtained 
straight line. This is well represented by our lin- 
ear fit which gives only a minor formal error. So, 
we think that despite the large formal errors of 
the astrometric fittings, the obtained f2i value, at 
least for the case in hand should be correct. 

We have to remark that the displacement 
of the radio source during our observing 
session was almost twice the value which 
was expected from the pure orbital mo- 
tion. Formally, of course, we were able to 
fit an astrometric orbit with a semi-major 
axis of aB = 0'.'0035, but the semi-major axis 
of the secondary's absolute orbit should be 
aB = 0'.'0019. Nevertheless, although our 
four-adjustable-parameter fit gave an un- 
reastically large value for the semi-major 
axis, and consequently, should be rejected, 
it gave the same value for ili as the three- 
parameter (fixed gb) fit. This also sug- 
gests, that despite the large formal errors, 
the value obtained for fli seems to be well- 
determined. This apparent large displace- 
ment or scatter is likely the consequence of 
the positional errors caused by short-term atmo- 
spheric phase fluctuations, and the variable struc- 
ture of the source during the flare. This would 
not be unprecedented. Large positional 
change was observed in the RS CVn sys- 
tem IM Peg during a flare by Lebach et al. 
(1999). These structural variations and the ori- 
gin of large radio flares in Algol could be studied 
with VLBI array configurations and observing fre- 
quencies providing (sub-)mas angular resolution. 
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4.2. Comparison to former VLBI mea- 
surements 

Before further discussion of the dynamical con- 
sequence of our result, we feel it necessary to com- 
ment on the well-known VLBI result of Lestrade 
et al. (1993). In our opinion it is without question 
that the excellent paper of Lestrade et al. (1993) 
is epoch-making in its significance, but, unfortu- 
nately, at the last step of their analysis they made 
some mistakes. As we cited earlier they obtained 
r^i = 52°. A careful look at their Fig. 3 clearly 
shows that this cannot be the correct result. One 
can see in that figure, that the coordinate differ- 
ence is larger in the declination direction than in 
the right ascension one. Consequently, the slope of 
the straight line fitted to their four points should 
be less than 45°, at least, when it is measured from 
north to east (i.e. from S to a). So, in our opin- 
ion, they obtained their value by measuring from 
east to north, and so their correct result should 
be r^i — 38° ± 5°. Furthermore, we found, 
that the exchange of the {S, a) coordinate 
pairs was not limited only to the determi- 
nation of rii, but it was applied in their all 
astro metric calculations. A less critical fur- 
ther consequence is that they obtained a 
reversed orbital revolution (compare their 
Fig. 4 with our Fig. [7]). However, the case of 
the correction for the orbital motion in the 
triple system is more problematic. Due to 
the aforementioned exchange of coordinate 
pairs, the direction of the orbital revolu- 
tion in the wide orbit is also reversed, and, 
consequently, the correction of the four ob- 
served coordinates for the orbital motion in 
the triple system is erroneous. 

Fig. [7] shows the corrected data points together 
with Lestrade et al.'s original solution. As one 
can see we obtain somewhat larger scatter in the 
data points. For these points we obtained Cli = 
45° ± 20° from the linear fit. (In this case 
we do not calculate an astro metric fit, as 
practically only two data points are known 
for the orbit. Remember, point one and 
two, as well as three and four belong almost 
to the same orbital phase, respectively.) 

We should note that in the case of the 
simple linear fits, the probable errors say 
nothing about the physical reliability of the 



results, or the accuracy of the measure- 
ments. They simply indicate the possibil- 
ity to fit one simple line for the four data 
points. To clarify this statement we have 
to keep in mind that the four points prac- 
tically belong to two orbital phases. Con- 
sequently, theoretically the first two points 
should practically coincide, and the same 
is also true for the third and fourth ones. 
Instead of this, one can see that the dis- 
tances of points one and two are w 0.9 mas 
and w 1.3 mas according to Lestrade's and 
our corrections, respectively, while for the 
other two points these values are « 1.0 mas 
and w 1.8 mas, respectively. These distances 
seem to be in good agreement with the 
statement of Mutel et al. (1998) about the 
radio source of Algol B, i.e. ,,The structure 
is double lobed with a separation of 1.6 ±0.2 
mas (1.4 times the K star diameter)" So 
this means that due to the extended, and 
presumably varying structure of the radio 
source, we cannot expect larger accuracy 
from the VLBI position measurement. Re- 
turning to the question of the probable er- 
rors, on using the Lestrade et al. (1993) 
original correction, the four data points 
then coincide almost in one straight line, so 
we can get a better linear fit than with our 
correction but as theoretically we should 
get only two points instead of four, this fact 
does not give any information about the re- 
liability of the two results. Furthermore, 
Mutel et al. (1998) conclude that the indi- 
vidual lobes are in the polar region, which 
is in better correspondance with the posi- 
tion of the radio source with respect to the 
orbit, in our ,,less accurate" solution (see 
again Fig. [7]). 

Taking into account the large scatter in 
the positions, this is in a very good agree- 
ment with the polarimetric measurements 
(f7i = 47° ± 7, Rudy 1979) and with our 
CHARA measurements (f^i = 48° ±3°) as 
well. 

4.3. Dynamics of the system 

In order to investigate the dynamical behaviour 
of Algol in the near past and future, we carried out 
numerical integration of the orbits for the triple 



8 



system. Detailed description of our code can be 
found in Borkovits et al. (2004). This code si- 
multaneously integrates the equations of the or- 
bital motions and the Eulerian equations of stel- 
lar rotation. The code also includes stellar dis- 
sipation, but the short time interval of the data 
allows that term to be ignored. Our input param- 
eters for Algol AB were almost identical with that 
of Table [T] with the exception of fii which was 
set to 48° in accordance with our CHAR A result. 
The orbital elements of the wide orbit were taken 
from Table[5](with two natural modifications, 
namely, instead of uab and luab, 0,2 and 012 
were used). As further input parameters, the 
^2, internal structure constants for the binary 
members were taken from the tables of Claret and 
Gimenez (1992) as /cf ^ = 0.0038, = 0.0011, 
fcf ^ = 0.0240, fcf ^ = 0.0087, respectively. 

Our numerical results between 1600 and 2100 
AD can be seen in Fig. [8] The variation of incli- 
nation between 7500 BC and 22 500 AD was also 
computed and can be seen in Fig. [9l Note that Al- 
gol AB does not show eclipses when the inclination 
is lower than 63° or higher than 117°. It shows 
partial eclipses if the inclination is between 63°- 
117° and moreover, it shows total eclipses when 
the inclination is between 87.3°~92.7. Therefore 
the last time when Algol was not an eclipsing bi- 
nary was before 161 AD and it showed partial 
eclipses between 161 AD and 1482 AD with in- 
creasing amplitude. Of course, at the beginning 
of this period, the eclipses featured a very small 
amplitude which later increased. By 1482 AD, the 
eclipses became total, and this was the case until 
1768 AD. The maximum length of the totality was 
about 0.5 hours around 1625 AD. It is wortnot- 
ing that in Algol the brighter star is the smaller 
one. Therefore the darker component could to- 
tally cover the brighter object causing large depth 
of minima of 2.8 magnitudes, so for a naked-eye 
observer it would almost disappear from the night 
for half an hour since its brightness during this half 
hour would be about 5.0 magnitudes. (The am- 
plitude nowadays is only about 1.3 magnitudes). 
Note, that during the totality the light of 
the wide, C component is the dominant. 
As one can see in these diagrams in the time of 
the discovery as a variable star (Montanari 1671), 
the inclination of the close pair was about 88°, 
making discovery easier. Our results also suggest 



that the light variation of Algol might have been 
known in the medieval Arabic and Chinese civil- 
isations (e.g. Wilk 1996). However, it should be 
emphasized that this time-data are rough approx- 
imations only. Since we could not determine the 
position of the node better than ±2°, and for ex- 
act calculations one needs an accuracy better by 
one order of magnitude, these numbers should 
be refined in the future. 

From about 1768 AD Algol show partial 
eclipses until approximately 3044 AD and the 
depth of the minima decreased in good agreement 
with the 20th century photometry measurements 
(Soderhjelm 1980). Considering the scientific era, 
one can see that our result suggests an inclination 
variation of Ai « — 1?6 in the last century which 
is in accordance with the statement of Soderhjelm 
(1980). 

5. Conclusions 

In this study we focused on the orientation of 
orbital planes in the hierarchical triple stellar sys- 
tem Algol. This is an important issue because 
the system has been showing a stable eclipse light 
curve for centuries: this could happen if the or- 
bital planes of the close pair and of the third 
body are either almost coplanar or perpendicu- 
lar to each other (Soderhjelm 1980, Borkovits et 
al. 2004). However, former estimations (Kiseleva 
et al. 1998 based on the results of Lestrade et al. 
1993) showed that the mutual inclination is 100° 
which would yield a fast inclination variation and 
consequently would result in the disappearance of 
the eclipses (Soderhjelm 1975, 1980; Borkovits et 
al. 2004). 

We found that joint use of optical and radio in- 
terferometry techniques in our project had great 
benefits. While in optical interferometry one can- 
not measure the visibility phase, there is a well 
understood a-priori source model (two stars orbit- 
ing each other). Because of this, even with the 
limited number of baselines available, we could fit 
well the value of the ascending node using visibility 
amplitudes only. In the radio regime, we are able 
to measure both visibility amplitude and phase, 
but we detect only the active corona of one of the 
stars. Unfortunately, this corona is highly vari- 
able and the Earth's atmosphere adds phase fluc- 
tuations that limit astrometric precision in short 
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measurements. However, with the combination 
of the VLBI measurements and the orbital 
phase information Eirising from the ecUpses 

one can resolve the ±180 phase ambiguity (once it 
is known which component emits in the radio). 

After careful analysis, we found the mutual in- 
clination angle of the orbital planes of the close 
and the wide pairs to be 95° ± 3°. Using this 
value as an initial value we integrated the equa- 
tion of motion of the system back to —7500 and 
forward to +22500. This helped to give support 
to the notion that medieval civilizations could ob- 
serve the big changes (up to 2.8 magnitudes) of 
Algol in the 17th century (Wilk 1996). The rate 
of inclination change of the close pair was found 
to be Ai w — l?6/century in the 20th century 
which shows only minor observable changes in the 
depth and shape of the minima in accordance with 
the photometric observations (Soderhjelm 1980). 
Therefore, the regular and precise observations of 
Algol's minima are recommended to further refine 
the geometrical configuration and to better under- 
stand the dynamics of triple stellar systems. 
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Table 1: Orbital elements, and astrophysical parameters of Algol A-B and AB-C determined by previous 
studies, ai, ei, ii and lui were taken from Kim (1989) while f^i is from Rudy (1979). The elements of the 
AB-C system were taken from Pan et al. (1993). The stellar quantities were also taken from Kim (1989). 



Quantity 


Notation 


A-B 


AB-C 


Time of periastron (HJD) 


T 


2445 739.0030^ 


2 446 931.4 


Period 


P 


2'?8673285 


680';'05 


Semi-major axis 


a 


0'.'0023'' 


0'.'09461 






14.1i?0 


582.9i?0'= 


Eccentricity 


e 




0.225 


Inclination 


i 


82?31 


83?98 


Argument of periastron 




_d 


310?29 


Longitude of the ascending node 


Q 


47° 


312?26 


Stellar Parameter 


Algol A 


Algol B 


Algol C 


Mass (Mq) 


3.8 


0.82 


1.8 


Radius (Rq) 


2.88 


3.54 


1.7 



"Time of primary minimum from Kim (1989). Note, if we set formally loi = 270° then this gives the time of periastron. 
''Kim (1989) gave the semi-major axis of the binary in solar radii (14.1i?Q). Using the HIPPARCOS parallax we transform it 
into arcseconds. 

■^Pan et al. (1993) gave the semi-major axis of the third body in arcseconds. Using the HIPPARCOS parallax we transform it 
into solar units. 

''Kim (1989) assumed a circular orbit. Hence eccentricity is zero and u) is not defined in the circular case. See also 
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Table 2: Log of observations and the observed normalized visibilities of Algol 





Time (\JT) 


n iinl 
U[IUJ 


V [lUJ 


B[m] 




V 


a 


{V) 


2006 Dec 2 




















W2-S2 


05:56:04 


54.709 


-168.117 


176, 


.794 


0, 


.763 


0, 


.031 


W2-S2 


06:01:56 


57.224 


-167.175 


176, 


.698 


0, 


.771 


0, 


.031 


W2-S2 


06:36:05 


71.063 


-160.881 


175, 


,877 


0, 


,749 


0, 


,029 


W2-S2 


06:59:09 


79.533 


-155.896 


175, 


.011 


0, 


.802 


0, 


.030 


W2-S2 


07:26:00 


88.366 


-149.424 


173, 


.597 


0, 


.738 


0, 


.030 


W2-S2 


07:48:43 


94.892 


-143.450 


171, 


,995 


0, 


,740 


0, 


,038 


E2-W2 


09:31:01 


60.207 


129.565 


142, 


.870 


0, 


.285 


0, 


.007 


E2-W2 


09:56:59 


44.772 


133.478 


140, 


.787 


0, 


.298 


0, 


.008 


E2-W2 


10:32:33 


22.735 


136.923 


138, 


,802 


0, 


,309 


0, 


,005 


E2-W2 


10:58:40 


6.173 


138.011 


138, 


.150 


0, 


.338 


0, 


.005 


E2-W2 


11:18:11 


-6.268 


138.009 


138, 


.151 


0, 


.419 


0, 


,003 


2006 Dec 3 




















E2-W2 


03:19:37 


130.402 


-0.093 


130.402 


0.616 


0.050 


E2-W2 


03:26:40 


132.347 


2.563 


132.372 


0.576 


0.050 


2006 Dec 4 




















E2-S2 


05:14:50 


-105.730 


-220.787 


244, 


,798 


0, 


,455 


0, 


,050 


E2-S2 


05:45:49 


-89.340 


-229.467 


246, 


.246 


0, 


.539 


0, 


.050 


W1-S2 


07:25:08 


196.777 


-142.326 


242, 


.854 


0, 


.601 


0, 


.020 


W1-S2 


07:35:15 


198.818 


-136.586 


241, 


,215 


0, 


,711 


0, 


,026 


W1-S2 


07:58:24 


202.018 


123.269 


236, 


.657 


0, 


.671 


0, 


.029 


W1-S2 


08:23:14 


203.147 


-108.827 


230, 


.461 


0, 


.623 


0, 


.021 


W1-S2 


08:46:36 


202.022 


-95.239 


223, 


,347 


0, 


,570 


0, 


,021 


W1-S2 


09:12:53 


198.237 


-80.137 


213, 


.823 


0, 


.560 


0, 


.041 


W1-S2 


09:33:28 


193.441 


-68.569 


205, 


.234 


0, 


.594 


0, 


.049 


W1-S2 


09:57:44 


185.776 


-55.361 


193, 


,850 


0, 


,608 


0, 


,035 


W1-S2 


10:19:34 


177.089 


-43.992 


182, 


.472 


0, 


.636 


0, 


.004 
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Table 3: Log of e-EVN observations and the observed total intensity peak positions of Algol. The complete 
data rows (following the horizontal lines) give the normal points formed from the each approx. 45-min. long 
observing scans, while in the case of the five-minutes averages only the astrometric angular coordinates are 
given. (See text for details.) (Algol position used for correlation (nominal coordinates): RA = 3^8™10!1315, 
DEC = +40°57'20'.'332 (J2000); [DEC = 40.955648, cos(DEC) = 0.7552172]; Fluxdcns: recovered VLBI 
total flux density during that timerange (mJy); R: modelfit centroid distance to phase centre (i.e. nominal 
position); 6: modelfit centroid position angle in degrees, measured from North to East; on the map North 
is top (Y axis) and East is to the left (X axis); X: measured nominal X coordinate in arcseconds ( 
X= 6RA ■ cos(DEC) ) X [arcsec] = R [mas]- sin 9/1000.0; SRA: X transformed to real RA coordinate 
difii'erence, in seconds; SRA [sec] = X [arcsec]/(15 • cos(DEC)); SDEC: measured - nominal Y coordinate in 
arcseconds (=^DEC); ^DEC [arcsec] = R[mas]- cose/1000.0) 



otarx \^ u ± J 


PnH /'TTT~\ 

LLiLLQ. ± J 


Fluxdsns. 






^ L J 




1 — OLJiijKj [ J 


iD:o/ 


I '7.1 Q 

I I :io 


4d.z mJy 


z4.oz 


iUo.4 




U.UUzll4 


-U.UUo / 1 


iD.O/ 


iD.OO 




fid 


1 no fi. 
iUZ.D 








iD.4i 


iD.4D 




OA OO 


1 HQ fi. 








io.4y 


ID. 04 




Z4.0D 


1 no o 
lUz.y 








iD.O / 


1 v. no 




OA QQ 


1 HQ Q 

iUo.y 








1 


1 7-1 n 




94 7Q 










17-1^ 
1 / .lo 


1 / . lo 




OA 

Z'i.OD 


lUo.U 








i ( :oo 


lo:io 


4o.4 mJy 


z4.oo 


iUo.y 


U.Uzooi 


n nnoi no 
U.UUziUz 


-u.uuosy 


1 V.QQ 
i / .OO 


1 7.QQ 
i / .OO 




OA Ad 
Z4.4D 










i 1 :4J. 


1 / :4o 




z4.4d 


Wo.6 








i ( :4y 


1 1 :J34 




zo.io 


iU4.1 








1 7-^7 


1 s-n9 




94 ^IS 
z^.ou 


1 04 ^ 








lo.UO 


lo.lU 




9/1 AO 
Z4.4Z 


lU.>.0 








1 o . 1 ■ ) 


1 ft- 1 ft 




Z i.. )Z 


l\) i.Z 








1 o .oo 

lo.oo 


1 n. 1 o 


oz.D nijy 


O /I /I o 

z4.4o 


1U4.4 


U.0^O<^ i 


0.002093 


-0.00609 


1 Q .QQ 


lo:oo 




O A PIT 

z4.o / 


1U4.D 








io.4i 


1o.4d 




OA 

Z4.04 


1 HQ O 

IUo.y 








1 Q./10 

io.4y 


io.04 




z4.0D 


1 f\A Q 

iU4.y 








18:57 


19:02 




24.31 


104.5 








19:05 


19:10 




24.42 


104.3 








19:13 


19:18 




24.45 


104.7 








19:33 


20:18 


25.0 mJy 


24.26 


105.2 


0.02341 


0.002067 


-0.00636 


19:33 


19:38 




24.26 


105.0 








19:41 


19:46 




24.12 


104.8 








19:49 


19:54 




24.43 


105.3 








19:57 


20:02 




24.33 


105.7 








20:05 


20:10 




24.32 


105.6 








20:13 


20:18 




24.07 


105.0 








20:33 


21:18 


16.5 mJy 


23.81 


105.5 


0.02294 


0.002025 


-0.00636 


20:33 


20:38 




23.78 


105.7 








20:41 


20:44 




24.38 


104.5 








21:01 


21:02 




24.02 


105.3 








21:05 


21:10 




23.71 


105.5 








21:13 


21:18 




23.74 


105.7 








21:33 


22:18 


12.6 mJy 


23.78 


106.2 


0.02284 


0.002016 


-0.00663 


21:33 


21:38 




23.62 


106.3 








21:41 


21:46 




24.01 


105.8 








21:49 


21:54 




23.90 


106.6 








21:57 


22:02 




23.85 


106.4 








22:05 


22:10 




23.44 


106.3 








22:13 


22:18 




23.85 


105.7 
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Table 3: Continuation of Table [3l 



Start (UT) End (UT) Fluxdens. R [mas] 9 X ["] 5RA [s] Y^^DEC ["] 



22:33 


22:18 


9.6 mJy 


23.51 


107.1 


0.02247 


0.001984 


-0.00691 


22:33 


22:38 




23.57 


106.7 








22:41 


22:46 




23.38 


106.8 








22:49 


22:54 




23.30 


107.6 








22:57 


23:02 




23.67 


107.5 








23:05 


23:10 




23.36 


107.6 








23:13 


23:18 




23.75 


106.9 








23:33 


00:18 


8.3 mJy 


22.92 


107.4 


0.02187 


0.001931 


-0.00685 


00:13 


00:18 




23.20 


107.2 








00:33 


01:18 


6.1 mJy 


22.90 


107.9 


0.02179 


0.001924 


-0.00704 


00:33 


00:38 




23.19 


107.7 








00:41 


00:46 




22.58 


107.6 








00:49 


00:54 




22.73 


108.9 








00:57 


01:02 




23.03 


109.0 








01:05 


01:10 




22.73 


108.5 








01:13 


01:18 




23.25 


105.7 
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Table 4: Results of the modeling of the CHARA observations.. N is the number of data points. 



Quantity 


Value 


Estimated uncertainty 


Surface brightness ratio in Kg band 


0.330 


±0.01 


Ascending node 


48° 


±2° 


Angular size of the semimajor-axis (mas) 


2.28 


±0.02 


XV{N-1) 


3.76 
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Table 5: Calciilatc'ci orbital dwiiwits of Algol B from EVN data. (Tlw^ first line at fii and belongs to the 
'fixed a' solution, while tln^ stx-ond one to the 'adjusted «' restilts. Errors art" la errors.) 



Quantity 


Notation 


Value 


Formal error 


Period 


Pi 


2?8673 


fixed 


Semi-major axis 


as 


0'/0019 


fixed 






c/ooasa 


0'.'00005 


Eccentricity 


ei 





fixed 


Inclination 


ii 


82?3 


fixed 


Argument of pcriastron 




Qa 


fixed 


Longitude of the ascending node 


Qi 


53° 


826° 






53° 


239° 


from linear fit 




52° 


3° 


Mean anomaly at to 


(Mo)i 


90° 


fixed 


Epoch 


to 


2 454084.360'' 






x' 


0.076623844 








0.025291438 




Quantities for the corrections 


Trigonometric parallax 


TT 


0'.'035 




Proper motion components 


Ma 


0!031cent-i 








-0'.'09cent-i 




Orl)ital elements of Algol AB in t 


(>rnar\- system'' 




Period 


P2 


679'?276349353 




Semi-major axis 


a-AB 


0'.'025738139 




Eccentricity 


62 


0.212719923 




Inclination 


«2 


84?014938082 




Argument of periastron 




132?558322883 




Longitude of the ascending node 


^2 


312?345789012 




Time of periastron 


T2 


2446 937.879685247 





"In case of circular orbit ui is undetermined. It was formally set to zero. This means that the mean anomaly (M) is measured 
from the ascending node. 
''Mid-eclipse moment of secondary minimum occurred during the EVN observation. 
'^Adopted from our recalculations of Pan et al. (1993) measurements. 
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"q^ Fig. 2. — The UV coverage of the Algol with 

L ff CHAR A. Each point represents one observational 

point. 
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Fig. 1. — The meaning of different angular orbital 
elements mentioned in this paper. Orbits are pro- 
jected into a sphere. Sli, 0.2 are the longitudes of 
the nodes of the close and the wide pair, respec- 
tively, A^i, N2 are the ascending nodes, ii and i2 
denotes the inclinations of the orbits measurable 
by an observer while im is the mutual inclination 
of the two orbital planes. Ji and J2 are the peri- 
centre points. 
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Fig. 3. — The result of the grid search for mini- 
mum value of ■ Only the result we obtained for 
ill is shown here. Note that the decimal loga- 
rithm of the /{N — 1) (N is the number of 
the data points) is shown for the sake of a better 
visualization. 
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Fig. 9. — The variation of the inchnation of Algol AB with time. Solid line represents the inclinitaion of the 
close pair observable from the Earth. Between the dashed lines the system shows partial eclipses, moreover 
between the dashed-dotted lines the system shows total eclipses for a short time. 
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Fig. 4. — The modelled versus observed visibility 
values and their errors. The solid line shows the 
1:1 line. 
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Fig. 5. — Radio total flux density variation of 

Algol from the Westerbork synthesis array data, 
taken during the VLBI observations. At the start 
of the observations there was a radio flare. 
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2 1.5 1 0.5 -0.5 -1 -1.5 -2 

Aa cosS (mas) 

Fig. 6. — The 5-minutes averaged data 
points with our astrometric fit of compo- 
nent B (see Table \5\ for the orbital ele- 
ments). The circled points were excluded 
from the fittings due to their large scatter. 
We labeled our first (16:45 heliocentric UT) 
and last (1:23 heliocentric UT) points. Fur- 
thermore, the 10*'S 20*'*, 30*'* and 40*'* points 
are also connected to their theoretical po- 
sitions on the orbit. The slim dashed line 
represents the line which connects the mid- 
eclipse points of the primary and the sec- 
ondary eclipses. 
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Fig. 7. — The astrometric orbit of Algol B 
for the four positions measured by Lestrade 
et al. (1993). Points 1 to 4 represents 
the points after Lestrade et al.'s original 
corrections, while points A to D denote 
our recalculated positions (see text for de- 
tails). (Note that the absolute coordinates 
of points 1 and A are not equal, but in this 
figure relative coordinates are used.) The 
observed points are connected by short lines 
with their corresponding positions along 
the astrometric orbits. The astrometric or- 
bits are also plotted: dashed lines repre- 
sents Lestrade et al.'s orbit while solid line 
corresponds to our recalculated orbit. The 
slim dashed lines represent the nodal line. 
(The ascending nodes are in the vicinity of 
points 2 and B, respectively). The arrows 
along the orbits show the direction of the 
orbital revolution. (Note, that the origi- 
nal solution of Lestrade et al. 1993 gives a 
reverse direction.) The arrows are located 
in the mid-eclipse point of the secondary 
minimum, where the mean anomly is equal 
to 90°. The inner of the two circles, cen- 
tered at the corresponding point along the 
recalculated orbit of position C represents 
the surface of the Algol B component, while 
the outer one illustrates schematically the 
separation of the two radio emitting lobes 
detected by Mutel et al. (1998). 
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Fig. 8. — The variation of the angular orbital ele- 
ments of Algol system between AD 1600 and 2100. 
The panels from top to bottom: {h) node of the 
binary measured in the invariable plane measured 
from the intersection of the invariable plane and 
the sky; (Qi) node of the close binary; (^2) same 
for the third component; (Zm) mutual inclination; 
(ii) inclination of the binary with respect to the 
invariable plane; (ii) inclination of the close bi- 
nary; (12) same for the tertiary. 
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ABSTRACT 

□ 

Algol is a triple stellar system consisting of a close semi-detached binary and a third object or- 
biting them. Due to the disputed spatial orientation of the close pair, the third body perturbation 

^Former address: Konkoly Observatory, H-1525 Budapest, P. O. Box 67. Hungary 
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of this pair is the subject of much research. In this study, we determine the spatial orientation 
of the close pair orbital plane using the CHARA Array, a six-element optical/IR interferometer 
located on Mount Wilson, and state-of-the-art e-EVN interferometric techniques. We find the 
longitude of the line of nodes for the close pair is Cldose = 48° ± 2° and the mutual inclination 
of the orbital planes of the close and the wide pairs is 95° ± 3°. This latter value differs by 5° 
from the formerly known 100° which would imply a very fast inclination variation of the system, 
not borne out by the photometric observations. We also investigated the dynamics of the system 
with numerical integration of the equations of motions using our result as an initial condition. 
We found large variation in the inclination of the close pair (its amplitude ~ 170°) with a period 
of about 20 millenia. This result is in good agreement with the photometrically observed change 
of the amplitude of Algol's primary minimum. 

Subject headings: stars: individual: Algol, stars: binaries 



1. Introduction 

There are about 1000 triple stellar systems 
known in the Galaxy, many of which consist of 
a close eclipsing pair and a distant third ob- 
ject orbiting around the close pair (Batten 1973, 
Tokovinin 1997). Algol is probably the most well- 
known such system. 

It consists of a semi-detached eclipsing binary 
with an orbital period of 2.87 days (B8V -|- K2IV) 
and an FlIV spectral type star revolves around 
them in every 680 days (discovered by radial ve- 
locity measurements, Curtiss 1908). Early inter- 
ferometric observations were unable to resolve the 
system (Merrill 1922), but the third component 
was succesfully observed by speckle interferometry 
(Gezari, Labeyrie and Stachnik 1972, Blazit et al. 
1977, McAlister 1977, 1979) and its orbit was pre- 
cisely determined by Bonneau (1979). This result 
was refined by using the Mark III optical stellar 
interferometer (Pan et al. 1993). 

In the radio regime Lestrade et al. (1993) de- 
tected positional displacement during the orbital 
revolution of the AB pair using the VLBI tech- 
nique, and identified the K-subgiant as the source 
of radio emission. The orbital elements of the close 
and the wide pairs determined from all these ob- 
servations are listed in Table 1. 

The light minima - mainly primary ones - of 
Algol were extensively observed in the last two 
centuries. There was only a very small change 
in the eclipse depth during this time. This led 
Soderhjelm (1975, 1980) to the theoretical con- 
clusion that the mutual inclination of the orbital 
planes of the close and the wide pair systems 
should not be larger than 11° and likely they are 
coplanar because both the shape and depth of the 



light minima should have noticeably changed oth- 
erwise. 

Since Pan et al. (1993) determined f^a = 312° 
for the third component (for the meaning of the 
various angles see Fig. ??), then the spatial orien- 
tation of the orbit has to be such that the longi- 
tude of the node of the the close pair in Algol 
should be about fli = 132° to have the coplanar 
configuration. Nevertheless, he noted that the po- 
larimetric result of Rudy (1979) gave f^i = 47°±7° 
for the node of the close orbit which would sug- 
gest a perpendicular rather than coplanar config- 
uration. 

This nearly perpendicular configuration was 
supported by other measurements: Lestrade et al. 
(1993) found Qi = 52° ± 5° degrees for the close 
pair in close agreement with the polarimetric mea- 
surements of Rudy (1979). If this value of f^i is 
correct, then the mutual inclination is about 100° 
(Kiseleva et al. 1998), and the two orbital planes 
are nearly but not exactly perpendicular to each 
other. This value for the mutual inclination of 
the system has been widely accepted since then. 
However, this mutual inclination value cannot be 
correct due to dynamical considerations. 

It is well-known that in a hierarchical triple 
stellar system, the orbital planes of the close and 
wide pairs are subject to precessional motion, in 
such a way that the normals of the orbital planes 
move on a conical surface around the normal of 
the invariable plane. (Here we omit the effect of 
stellar rotation which is insignificant in an ordi- 
nary triple system.) In the case of a hierachical 
triple system where the invariable plane almost 
coincides with the wider orbital plane, the preces- 
sion cone angle is close to the mutual inclination. 
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Consequently, in the case of the present system, we 
would get an almost 160°-amplitudc variation in 
the observable inclination during the approximate 
period given by Eq. (27) of Soderhjelm (1975) (as- 
suming the present approximation is valid, as the 
mutual inclination tends to 90°, the precession pe- 
riod tends to infinity). Here we refer to Fig. 4 of 
Borkovits et al. (2004) which clearly shows that in 
the case of the above mentioned configuration the 
observable inclination of the close binary would 
have changed approximately 3° in the last cen- 
tury which evidently contradicts the observations. 
In this case the eclipses would disappear within a 
few centuries. This was already observed in some 
eclipsing binaries, like in SS Lacertae or V907 Sco 
(see e.g. Eggleton & Kiseleva-Eggleton 2001) but 
not in Algol. 

In summary, the polarimetric and the inter- 
ferometric observations contradict the coplanar- 
configuration, but since the mutual inclination 
is far from the exact perpendicularity - they are 
not in agreement with the observed tiny change 
in the minima depth cither. A closer approximate 
perpendicularity of the two orbital planes would 
mean that the period of orbital precession becomes 
so large that the inclination variation (and con- 
sequently the depth variation of the minimum) 
of the close pair remains unobservable for a long 
time, consistent with the observations. 

The aim of this study was to constrain the mu- 
tual inclination of the system better, requiring the 
measurement of the longitude of the node for the 
close pair. The other orbital elements are well- 
known from spectroscopic or photometric data, 
but there is a controversy in the value of fli. 
Because the expected apparent size of Algol- AB 
semi-major axis is of the order of 2 milliarcseconds 
(mas), we carried out optical and radio interfer- 
ometry measurements. As we will show, optical 
and radio interferometry are complementary tech- 
niques. By joining these it is possible to resolve 
the ambiguity in the geometry of the system, and 
better assess the accuracy of our measurements. 

2. Observations and data reduction 
2.1. CHARA Observations 

The CHARA Array is an optical/ncar-IR in- 
terferometer array consisting of six Im telescopes. 
The array is described in detail in ten Brumme- 



laaretal. (2005). A detailed overview and further 

references about the obscrvablcs and the theory 
of optical interferometry can be found in Haniff 
(2007). 

Wc observed Algol on three nights (2, 3 and 4 
December, 2006) in Kg band. Much of the sec- 
ond night was lost due to high winds and dusty 
conditions. 

Iota Persei and Theta Persei served as calibra- 
tion stars. The observations of the target and the 
two calibrators were organized into a sequence and 
the measurements on calibrators generally brack- 
eted the target observations. We have 12 data 
points of Iota Persei, 12 data points of Theta Per- 
sei and 23 data points of Algol itself. 

Each data point was calculated from a number 
of scans. During a scan the intensity variations 
function of the path delay were measured. 
About 300 scans were collected within five min- 
utes for one data point, the first 22 of them were 
obtained on the targets (Algol or one of the cali- 
brator stars), then 21 scans were done for measur- 
ing the background while the shutter was closed, 
then more than 200 other scans were obtained on 
the targets again and finally 67 further scans for 
measuring the background, again. 

To reduce the data, wc used the recipe of McAl- 
ister (2002). This consisted of the following steps. 
First a low-pass filter was applied to remove the 
atmospheric noise. Then the bias was subtracted 
and the scans were normalized to unity. As a 
next step, the scans measured by the two chan- 
nels were subtracted from each other (for details, 
see ten Brummelaar et al. 2005 and McAlister 
2002). This was further processed by applying 
a high-frequency filter to reduce noise. Discrete 
Fourier-transforms of the scans were calculated 
and a template was computed. For this new tem- 
plate, we used the full amplitude for the frequen- 
cies ±25/cycle around maximum frequency and 
20% of the amplitude for the other frequencies. 
From this wc could calculate the maximum devi- 
ation of the template from zero which yielded an 
estimation of the visibility value. Removing out- 
lier values, we averaged the remaining ones which 
yielded the uncalibrated visibility value of this 
point. The error was estimated as the standard 
deviations of the visibility values of the more than 
200 scans of the point. 
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The measured visibilities of the calibrators were 
linearly interpolated for the times of the Algol ob- 
servations. A comparison of the true and mea- 
sured visibilities of the calibrators yielded a factor 
which converted the measured target visibilities to 
true ones. The true visibilities of the two calibra- 
tors were estimated as follows. 

The uniform disc (UD) angular diameter of Iota 
Persei is 1.21 ± 0.06 mas according to 'Catalog of 
High Angular Resolution Measurement' (Richichi 
et al. 2005). The true diameter of the other cali- 
brator star Theta Persei is 1.201 ± O.OlSi?© which 
was determined by the fit of its spectrum (Valenti 
and Fisher 2005). Since its HIPPARCOS parallax 
is known one can easily calculate that its true UD 
angular diameter is 0.995 ± 0.01 mas. 

The limb-darkened (LD) angular diameter of 
a stars is larger than the UD diameter. There 
is a simple relationship between them (Hanbury 
Brown et al. 1974): 

Old ^ / 1 - ma/3 , , 

OuD ^ 1- 7ux/15 ^ ' 

The limb- darkening coefficient were taken from 
the tables of van Hamme (1993). These coeffi- 
cients are a function of surface gravity and effec- 
tive temperature which were estimated from the 
known spectral type of the calibrators. Then Eq. 
(1) yielded the corrections which increase the UD 
angular diameters by a few percent only. These 
corrected values were used to calibrate the visibil- 
ities. 

The telescope combinations, epoch of observa- 
tions, baselines, uncalibrated visibilities and their 
errors can be found in Table 2. The UV-coverage 
can be seen in Fig. ?? 

2.2. e-VLBI Observations 

We observed Algol with a subset of the Eu- 
ropean VLBI Network (EVN) on 14-15 Decem- 
ber 2006 at 5 GHz for 9 hours. These obser- 
vations were carried out using the e-VLBI tech- 
nique, where the telescopes stream the data to 
the central data processor (JIVE, Dwingeloo, the 
Netherlands) instead of recording. Using e-VLBI 
for the observations was not fundamental for our 
measurements, but we took the opportunity of an 
e-VLBI test close in time to the CHARA obser- 
vations, outside the normal EVN observing ses- 



sion. The participating telescopes were Cam- 
bridge and Jodrell Bank (UK). Medicina (Italy), 
Onsala (Sweden), Toruh (Poland) and the West- 
erbork phased array (the Netherlands). The data 
rate per telescope was 256 Mbps, which resulted 
in 4x8 MHz subbands in both LCP and RCP po- 
larizations using 2-bit sampling. The correlation 
averaging time was 2 seconds, and we used 32 de- 
lay steps (lags). Initial clock searching was car- 
ried out before the experiment using the fringe- 
finder source 3C345. Algol was phase-referenced 
(Beasley & Conway 1995) to 0309-^411 in 3-5-3 
minute cycles. Additional scans were scheduled on 
3C84 for real-time fringe monitoring, and for D- 
term calibration. We used 3C138 to cahbrate the 
Westerbork synthesis array amplitudes and polar- 
ization. 

Post-processing was done using the US Na- 
tional Radio Astronomy Observatory (NRAO) 
AIPS package (Diamond 1995). The amplitudes 
were calibrated using the known antenna gain- 
curves and the measured system temperatures. 
The data were fringe-fitted, bandpass calibrated, 
and then polarization calibrated. We corrected 
for the polarization leakage D-terms and fringe- 
fitted the cross-hand data, after which the data 
were averaged in frequency in each subband. Be- 
sides the standard procedure, we used the WSRT 
synthesis array measurements on 0309+411 to 
obtain a more accurate VLBI flux scale. The 
phase-reference source showed a low level of cir- 
cular polarization (fractional CP ^ 0.28%). The 
left and right-handed VLBI gains were separately 
adjusted in accordance with the WSRT mea- 
surement. Imaging was carried out in Difmap 
(Shepherd et al. 1994). The snap-shot images 
(from about 46 minutes data each) were made by 
Fourier-transforming the observed visibilities, no 
self-calibration was applied. The position of Al- 
gol was thus determined within each snap-shot by 
directly model-fitting the uw-data in Difmap. 

3. Data analysis and results 
3.1. Analysis of CHARA data 

According to the van Cittert-Zernike theorem, 
the amplitude of the visibility is the normalized 
Fourier-transform of the intensity distribution (for 
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a comprehensive explanation, see Haniff 2007): 



V{u,v,t) 



J J I{x,y,t)cos(2Tr{u{t)x + v{t)y)/ \)dxdy 
J J y, t)dxdy 



(2) 

In this equation V is the true visibihty at the {u, 
v) spatial frequencies, {x, y) are the corresponding 
sky coordinates, t is the time, / is the intensity at 
the {x,y) sky point and finally A = 2.133 fjtm. is 
the central wavelength of our observations (corre- 
sponding to the Ks band). Note that the inten- 
sity distribution on the sky normally changes very 
slowly with time, but we had to include the time 
dependence in Eq. (2) because of the rapid orbital 
revolution of the AB pair. Because of the power 
of CHARA, it is not enough to model the system 
as the combination of two point-like sources but 
we need to combine the pictures of two extended 
sources. 

To determine this intensity distribution wc de- 
veloped a model which was very close to the one 
of Wilson and Devinney (1971) which is based on 
the Roche-model (Kopal 1978) but it was imple- 
mented in IDL to restore the surface intensities 
into a matrix and to calculate the sky-projected 
picture of the system. 

Since we have only about two dozen visibility 
measurements, we wanted to limit the number of 
free parameters, choosing just three: the angu- 
lar size of the semi-major axis (in mas), the Ks 
surface brightness ratio, and the angle ili. Every 
other parameters were fixed according to the val- 
ues given in Wilson et al. (1972) or in Kim (1989). 

The model outlined above was used to fit the 
data and a grid search was carried out on the free 
parameters: 

a) The surface brigthness ratio was stepped 
from to 1.4 with a step size of 0.05; 

b) rii was stepped from 0° to 180° with a step- 
size of 5° (note that the visibility amplitude 
does not change if we rotate the image by 
180°, which leaves a 180° ambiguity in the 
node of the periastron; this ambiguity in Oi 
can be resolved with VLBI); and 

c) the angular size of the semi-major axis was 
stepped from 1.80 mas to 3.60 mas with a 
step-size of 0.005 mas (note that the ex- 
pected size was 2.50 mas). 



Nearly 353 000 models were calculated on this 
grid, and the minimum was found. Around the 
minimum, a new search was carried out with a 
finer grid and about 35 000 new models were com- 
puted again. Around the minimum, a polynomial 
fit yielded the final values and errors. The results 
are shown in Figs. ?? and ??. 

The best solution we found is reported in Ta- 
ble 3. 

3.2. Analysis of e-VLBI data 

The EVN measurements were carried out be- 
tween 16:39 and 1:19 U.T. on 14/15 December 
2006. As the orbital period of the eclipsing bi- 
nary is P = 2.8673 days, the covered orbital arc 
was near 13 %. During this period a secondary 
minimum occurred which was observed simulta- 
neously photometrically by the 50-cm telescope at 
Piszkesteto Station of the Konkoly Observatory, 
Hungary (see Biro et al., 2007). One of our goals 
was to observe possible partial occultation of the 
radio source by the primary component and mea- 
sure the change in the circular polarization proper- 
ties of the source accordingly. Because Algol flared 
during the observations, this goal could not be ful- 
filled and will not be discussed further here. 

In contrast to the measurements of Lestradc ct 
al. (1993) who observed the Algol on four differ- 
ent nights in near quadrature phases (i.e. around 
<j} = 0.25 and (p ~ 0.75), our observations cov- 
ered a small part of the orbit around (j) ~ 0.5, 
when the arc projected onto the plane of the sky 
is the largest. (Naturally the case is the same 
at ~ 0). The advantage of this approach is 
that the motion of the target can be detected in 
a few hours, and the measured positions are only 
slightly affected by the orbital motion of the AB-C 
pair, unlike in the case when the data are taken 
at different epochs. Because of the short observ- 
ing time interval, the orbit of the AB pair itself 
is not well constrained. There are two limitations 
that must be mentioned here, 1) short-term tropo- 
spheric and ionospheric phase fluctuations which 
cannot be modelled well and limit the astromet- 
ric accuracy in short VLBI-mcasurements, and 2) 
the variable structure of the radio source (Mutel 
et al. 1998) - the radio emission is not coming 
from the K-subgiant, but likely from its active po- 
lar region. Although Algol was unresolved with 
our array configuration, the source centroid po- 
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sition could have changed appreciably because of 
the bright flare during the run. These are the fac- 
tors that must be taken into account in the inter- 
pretation of the final result. 

We calculated astrometric orbits as well 
simple linear fit on two different sets of observing 
data. First, hourly normal points were formed. 
Then we also calculated the orbit using 6 minutes 
averages. The latter showed that some points with 
extremely large scatter produce false result in the 
hourly normal points. Removing the outliers, wc 
calculated our final solution from the 6-minute- 
average data. Despite the shortness of our observ- 
ing session (less than 9 hoius) the positional data 
were corrected for the annual parallax, proper mo- 
tion and the revolution around the centre of mass 
of the triple system. These corrections resulted in 
approximately 1° difference in the node position. 
For this latter correction we recalculated the third 
body orbit by the same code which was used for 
the binary orbit determination from our e-VLBI 
data. We used both the data sets of Bonneau 
(1979) and Pan et al. (1993). From these (very 
similar results) we applied the orbital elements 
obtained from Pan's data (see Table ??) for the 
wide orbit correction. For these calculations, we 
used our own differential correction code based on 
a Levenberg-Marquardt algorithm, tested against 
the data and results of Eichhorn and Xu (1990), 
and found it to be in excellent agreement. The 
calculation consists of nine free parameters: the 
position of the centre of mass (Xo, Ya), the orbital 
period (P), and the six usual orbital elements (a, 
e, i, ui, Oi, and Mq). Because of the circular orbit, 
instead of the argument of periastron (w), and the 
mean anomaly (M), their sum, i.e., the true lon- 
gitude (the distance from the ascending node in 
case of circular orbits) should be used. This was 
done in such a way, that u) was formally consid- 
ered as zero, while Mq = 90° was set for the mid- 
(secondary) eclipse moment, to = 2 454 084.360. 
The period and the inclination was acquired from 
Kim (1989), while the semi-major axis of the sec- 
ondary's orbit around the centre of mass of the 
binary was calculated from Kim's data. First we 
adjusted three parameters {Xq, Yq, Cli), leaving 
the other six fixed. Finally, we make adjustments 
to the semi-major axis. 

Using Gnuplot, we also performed a simple lin- 
ear fit to the data from the knowledge that in the 



vicinity oi (p = 0.5 the motion can be approxi- 
mated well with a line whose slope gives tanfJi. 

Our solutions can be seen in Table ??, as well 
as in Fig. ??. 

4. Discussion 

4.1. Discussion of the results 

First, we concentrate on the CHARA results. 
The true size of the semi-major axis of Algol is 
14.1 Rq (Kim 1989) while its angular size was 
measured by us to be 2.28 ± 0.02 mas (see Ta- 
ble 3). By dividing these two numbers one can 
find that the distance to Algol is 28.6 ± 0.3 par- 
sec. The HIPPARCOS parallax yielded 28.46t^-j5 
parsec. The agreement is excellent. 

The determined surface brightness ratio (0.33) 
can be converted to a luminosity ratio by multi- 
plying it with the ratio of surface area of the two 
stars, which is known from the light curve solu- 
tion (Wilson et al. 1972; Kim 1989). The resulted 
luminosity ratio is 0.43. This value is very close 
to the photometrically estimated 0.44 (Murad and 
Budding 1984). 

According to the CHARA results, the angle of 
the node is f2i = 48° ± 2° with an ambiguity of 
180 degrees. Because the determined distance and 
luminosity ratios agree very well with earlier mea- 
surements obtained with other methods, we are 
confident in our results. With the VLBI measure- 
ments (see below) we resolve the ±180° ambiguity 
and conclude that Oi =48° ±2°. This is in ex- 
cellent agreement with the value determined from 
polarimetric measurements (fii = 47° ± 7°, Rudy 
1979), indicating that polarimetry is an efficient 
tool to determine the spatial orientation of the or- 
bits. 

At this point we can determine the mutual in- 
clination im with the following formula: 

cosim = cosii cos 12 -l- sinzi sin 12 cos(f2i — SI2) 

(3) 

The result is ifn — 95° it 3° (the uncertainty re- 
flects to the uncertainties not only in fli but ones 
in the other angular elements), confirming the con- 
clusion of Lestrade et al. (1993) that the two 
orbital planes are nearly perpendicular to each 
other. This value is, however, closer to the exact 
perpendicularity than the 100° given in Kiseleva et 
al. (1998) which was based on the measurements 
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of Lestrade et al. (1993). Nevertheless, the exact 
perpendicularity is within the three sigma range. 

Regarding the e-VLBl measurements (see Ta- 
ble 4), we got Oi = 53° ± 239° for the node from 
both astromctric fits, and fli ~ 52° ± 3° from tlic 
simple linear one. These values are close to the 
previously obtained ones. However, we note, that 
in the case of the astrometric orbit fitting, the for- 
mal errors are extremely large. This naturally re- 
flects the fact, that our measurements cover only 
a very short fraction of the orbit, and especially in 
that phase, where the expected astrometric orbit 
is almost a straight line. Consequently, without 
any a priori information, the orbit would be com- 
pletely undetermined. However, in this particular 
case, the longitude of the node itself is very well 
determined during this phase, as this is nothing 
other than the slope of the obtained straight line. 
This is well represented by our linear fit which 
gives only a minor formal error. So, we think that 
despite the large formal error of the astrometric 
fittings, the obtained Oi value, at least from that 
point of view, should be correct. 

Now we turn our attention to the velocity of 
the radio source. As can be easily calculated from 
the data listed in Table ??, the semi-major axis of 
the secondary around the barycenter of the close 
binary should be 0'.'0019, while the observed arc 
belongs to an orbit with a = 0'.'0035. This may 
simply be the consequence of the positional er- 
rors caused by short-term atmospheric phase fiuc- 
tuations, or perhaps the variable structure of the 
source during the large flare. These structural 
variations and the origin of large radio flares in 
Algol could be studied with VLBI array conflgu- 
rations and observing frequencies providing (sub- 
)mas angular resolution. 

4.2. Comparison to former VLBI mea- 
surements 

Before further discussion on the dynamical con- 
sequence of our result, we feel it necessary to com- 
ment the well-known VLBI result of Lestrade et al. 
(1993). In our opinion it is out of question that 
the excellent paper of Lestrade et al. (1993) is an 
epoch-making one in its significance, but, unfortu- 
nately, at the last step of their analysis they made 
some mistakes necessary to correct. As we cited 
earlier they obtained f2i = 52°. A careful look on 
their Fig. 3 clearly shows that this cannot be the 



correct result. One can see in that figure, that the 

coordinate difference is larger in the declination 
direction than in the right ascension one. Conse- 
quently, the slope of the straight line fitted to their 
four points should be less than 45°, at least, when 
it is measured from north to east (i.e. from 6 to 
a). So, in our opinion, they got their value mea- 
suring it from east to north, and so their correct 
result should be f^i = 38° ± 5°. 

In that case, of course, the question of the dis- 
crepancy between the earlier polarimetric, as well 
as our CHARA data, and their VLBI result nat- 
urally arises. In our opinion, the answer can be 
found in the data position correction. As their 
four data points were collected in a more-than- 
seven-day-long interval, the correction for the an- 
nual parallax and the orbital motion around the 
centre of mass of the ternary system is critical. We 
tried to recalculate the corrected values of their 
measurements departing from their original data, 
and we found that in their work there is a 180° dis- 
crepancy in the argument of periastron, and, con- 
sequently, in the position of the binary in its wide 
orbit. The corrected data, after removing this dis- 
crepancy can be seen in Fig. ??. In that case we 
got somewhat larger scatter in the data points. 
For these points we obtained Qi = 45° ± 20° and 
fii = 56° ± 429° from the linear and the astro- 
metric fits, respectively. Taking into account the 
large scatter in the positions, these results are in 
accordance with the other measurements. Since a 
simple linear fit has much smaller uncertainty one 
can accept this value (fii = 45° ±20°) as a correct 
result of Lestrade et al. (1993). Then this also is in 
a very good agreement with the polarimetric mea- 
surements {fli = 47° ± 7, Rudy 1979) and with 
our CHARA measurements (ili = 48° ± 3°). 

4.3. Dynamics of the system 

In order to investigate the dynamical behaviour 
of Algol in the near past and future, we carried out 
numerical integration of the orbits for the triple 
system. The detailed description of our code can 
be found in Borkovits et al. (2004). This code 
simultaneously integrates the equations of the or- 
bital motions and the Eulerian equations of stel- 
lar rotation. The code also includes stellar dis- 
sipation, but the short time interval of the data 
allows that term to be ignored. Our input pa- 
rameters for Algol AB were almost identical with 
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that of Table ?? with the exception of fii which 
was set to 48° in accordance with our CHARA 
result. The orbital elements of the wide orbit 
were taken from Table ?? (with a natural 180° 
shift in Q2)- As further input parameters, the 
k2, internal structure constants for the binary 
members were taken from the tables of Claret and 
Gimenez (1992) as /cf ^ = 0.0038, = 0.0011, 
ki^^ = 0.0240, = 0.0087, respectively. 

Our numerical results between 1600 and 2100 
AD can be seen in Fig. ??. The variation of incli- 
nation between -7500 BC and 22 500 AD was also 
computed and can be seen in Fig. ??. Note that 
Algol AB does not show eclipses when the incli- 
nation is lower than 63° or higher than 117°. It 
shows partial eclipses if the inclination is between 
63°-117° and moreover, it shows total eclipses 
when the inclination 87.3°-92.7. Therefore the 
last time when Algol was not an eclipsing binary 
was before 24 AD and it showed partial eclipses 
between 24 AD and 1345 AD with increasing am- 
plitude. Of course, at the beginning of this period, 
the eclipses were ones with a very small ampli- 
tude, but later increased in amplitude. By 1345 
AD, the eclipses became total, and this was the 
case until 1631 AD. The maximum length of the 
totality was about 0.5 hours around 1488 AD. It 
is worthy to note that in Algol the brighter star 
is the smaller one. Therefore the darker compo- 
nent could totally cover the brighter object caus- 
ing as large depth of minima as 2.8 magnitudes 
so for a naked-eye observer it almost disappeared 
from the night for half hour since its brightness 
during this half hour was about 5.0 magnitudes. 
(The amplitude nowadays is only about 1.3 mag- 
nitudes). As one can see in these diagrams in the 
time of the discovery as a variable star (Montanari 
1671), the inclination of the close pair was about 
88° which yielded an amplitude of 1.6 magnitudes, 
making discover easier. It also is evident that our 
results arc also in accordance with the fact, that 
the light variation of Algol might have been known 
in the medieval Arabic and Chinese civilisations 
as well (e.g. Wilk 1996). However, it should be 
emphasized that this time-data are rough approx- 
imations only. Since we could not determine the 
position of the node better than 2° and for en ex- 
act calculations one needs a better accuracy with 
one order of magnitude, these numbers should be 
refined in the future. 



Prom about 1631 AD it showed partial eclipses 

until approximately 3044 AD and the depth 
of the minima decreased in good agreement 
with the 20th century photometry measurements 
(Sodcrhjelm 1980). Considering the scientific era, 
one can see, that our result suggests an inclination 
variation of Ai » — 1?6 in the last century which 
is in accordance with the statement of Soderhjelm 
(1980). 

5. Conclusions 

In this study we focused on the orientation of 
orbital planes in the hierarchical triple stellar sys- 
tem Algol. This is an important issue because 
the system has been showing a stable eclipse light 
curve for centuries: this could happen if the or- 
bital planes of the close pair and of the third 
body arc either almost coplanar or perpendicu- 
lar to each other (Soderhjelm 1980, Borkovits et 
al. 2004). However, former estimations (Kiseleva 
et al. 1998 based on the results of Lcstrade et al. 
1993) showed that the mutual inclination is 100° 
which would yield a fast inclination variation and 
consequently would result in the disappearance of 
the eclipses (Soderhjelm 1975, 1980; Borkovits et 
al. 2004). 

We find that joint use of optical and radio in- 
terferometry techniques in our project had great 
benefits. While in optical interferometry one can- 
not measure the visibility phase, there is a well un- 
derstood a-priori source model (two stars orbiting 
each other). Because of this, even with the lim- 
ited number of baselines available, we could well 
fit the value of the ascending node using visibility 
amplitudes only. In the radio regime, we are able 
to measure both visibility amplitude and phase, 
but we detect only the active corona of one of the 
stars. Unfortunately, this corona is highly vari- 
able and the Earth's atmosphere adds phase fluc- 
tuations that limit astrometric precision in short 
measurements. However, with VLBI one can re- 
solve the ±180 phase ambiguity (once it is known 
which component emits in the radio). 

After careful analysis, we found the mutual in- 
clination angle of the orbital planes of the close 
and the wide pairs to be 95° ± 3° at epoch 
2007.92. Using this value as an initial value we in- 
tegrated the equation of motion of the system back 
to -7500 and forward to -1-22500. This helped to 
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give support that medieval civilizations could ob- 
serve the big changes (up to 2.8 magnitudes) of 
Algol in the 15th century (Wilk 1996). The rate 
of inclination change of the close pair was found 
to be Ai « — l?6/century in the 20th century 
which shows only minor observable changes in the 
depth and shape of the minima in accordance with 
the photometric observations (Soderhjelm 1980). 
Therefore, the regular and precise observations of 
Algol's minima are recommended to further refine 
the geometrical configuration and to better under- 
stand the dynamics of triple stellar systems. 
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Table 1: Orbital elements, and astrophysical pa- 
rameters of Algol A-B and AB-C determined by 
previous studies, a, e, i and oj were taken from 
Kim (1989) while Qi is from Rudy (1979). The el- 
ements of the AB-C system were taken from Pan ct 
al. (1993). The stellar quantities were also taken 
from Kim (1989). 



A-B 


AB-C 


).0030^ 


2 446 931.4 


73285 


680^05 


|'.'0023^ 


0'.'09461 


LIRq 


582.9i?o'= 




0.225 


82?31 


83?98 


_d 


310?29 


47° 


312?26 


Igol B 


Algol C 


0.82 


1.8 


3.54 


1.7 



"Time of primary minimum from Kim (1989). Note, if we set formally ui = 270° then this gives the time of periastron. 
''Kim (1989) gave the semi-major axis of the binary in solar radii (14.I0). Using the HIPPARCOS parallax we transform it into 
arcseconds. 

°Pan et al. (1993) gave the semi-major axis of the third body in arcseconds. Using the HIPPARCOS parallax we transform it 

into solar units. 

''Kim (1989) assumed a circular orbit. Hence eccentricity is zero and u) is not defined in circular case. See also a. 
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Table 2: Log of observations and the observed nor- 
malized visibilities of Algol 



Telescopes 


Time (UT) 


u[m] 


v[m] 


B[m] 




V 


a 


iV) 


2006 Dec 2 


















W2-S2 


05:56:04 


54.709 


-168.117 


176.794 


0, 


.763 


0, 


.031 


W2-S2 


06:01:56 


57.224 


-167.175 


176.698 


0, 


,771 


0. 


.031 


W2-S2 


06:36:05 


71.063 


-160.881 


175.877 


0, 


.749 





.029 


W2-S2 


06:59:09 


79.533 


-155.896 


175.011 


0, 


.802 


0, 


.030 


W2-S2 


07:26:00 


88.366 


-149.424 


173.597 


0, 


.738 


0, 


.030 


W2-S2 


07:48:43 


94.892 


-143.450 


171.995 


0, 


.740 


0, 


.038 


E2-W2 


09:31:01 


60.207 


129.565 


142.870 


0, 


.285 


0, 


.007 


E2-W2 


09:56:59 


44.772 


133.478 


140.787 


0, 


.298 


0, 


.008 


E2-W2 


10:32:33 


22.735 


136.923 


138.802 


0, 


.309 


0, 


.005 


E2-W2 


10:58:40 


6.173 


138.011 


138.150 


0, 


.338 


0, 


.005 


E2-W2 


11:18:11 


-6.268 


138.009 


138.151 


0, 


.419 


0. 


.003 


2006 Dec 3 


















E2-W2 


03:19:37 


130.402 


-0.093 


130.402 


0.616 


0.050 


E2-W2 


03:26:40 


132.347 


2.563 


132.372 


0.576 


0.050 


2006 Dec 4 


















E2-S2 


05:14:50 


-105.730 


-220.787 


244.798 


0, 


.455 


0, 


.050 


E2-S2 


05:45:49 


-89.340 


-229.467 


246.246 


0, 


.539 


0, 


.050 


W1-S2 


07:25:08 


196.777 


-142.326 


242.854 


0, 


.601 


0. 


.020 


W1-S2 


07:35:15 


198.818 


-136.586 


241.215 


0, 


.711 


0, 


.026 


W1-S2 


07:58:24 


202.018 


123.269 


236.657 


0, 


.671 


0, 


.029 


W1-S2 


08:23:14 


203.147 


-108.827 


230.461 


0, 


.623 


0, 


.021 


W1-S2 


08:46:36 


202.022 


-95.239 


223.347 


0, 


.570 


0, 


.021 


W1-S2 


09:12:53 


198.237 


-80.137 


213.823 


0, 


.560 


0, 


.041 


W1-S2 


09:33:28 


193.441 


-68.569 


205.234 


0, 


.594 


0, 


.049 


W1-S2 


09:57:44 


185.776 


-55.361 


193.850 


0, 


.608 


0, 


.035 


W1-S2 


10:19:34 


177.089 


-43.992 


182.472 


0, 


.636 


0, 


.004 
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Table 3: Results of the modeling of the CHARA 
observations. 



Quantity 


Value 


Estimated uncertainty 


Surface brightness ratio in band 


0.330 


±0.01 


Node of periastron 


48° 


±2° 


Angular size of the semimajor-axis (mas) 


2.28 


±0.02 




3.76 
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Tabic 4: Calculated orbital elements of Algol B 
from EVN data. (The first row at il and be- 
longs to the 'fixed a' solution, while the second 
one to the 'adjusted o' results.) 



Quantity 


Notation 


Value 


Formal error 


Period 


P 


2'?8673 


fixcxl 


Semi-major axis 


a2 


0'.'0019 


fixed 






0'.'00353 


0'.' 00005 


Eccentricity 


e 





fixed 


Inclination 


i 


82?3 


fixed 


Argument of periastron 




Qa 


fixed 


Longitude of the ascending node 


n 


53° 


826° 






53° 


239° 


from linear fit 




52° 


3° 


Mean anomaly at to 


Mo 


90° 


fixed 


Epoch 


to 


2 454 084.360'^ 






x' 


0.076623844 








0.025291438 




Quantities for the corrections 


Trigonometric parallax 


IT 


C.'OSS 




Proper motion components 




W031cent-i 






I^S 


-0'.'09cent-i 




Orbital elements of Algol AB in ternary system*^ 


Period 


P' 


679^276349353 




Semi-major axis 


a'i2 


0'.'025738139 




Eccentricity 


e' 


0.212719923 




Inclination 


i' 


84?014938082 




Argument of periastron 


'^12 


132?558322883 




Longitude of the ascending node 


fl' 


312? 345789012 




Time of periastron 


T' 2 446 937.879685247 





"In case of circular orliil uj is undetermined. It was formally set to zero. It means, that the mean anomaly (M) is measured 
from the ascending node. 

''Mid-eclipse moment of secondary minimum occurred during the EVN observation. 
^Adopted from our recalculations of Pan et al. (1993) measurements. 



14 



300 



To Observer 




200 



100 



^ 



-100 



-200 



-300 



East 



Fig. 1. — The meaning of different angular orbital 
elements mentioned in this paper. Orbits are pro- 
jected into a sphere, fii, il.2 are the longitudes of 
the nodes of the close and the wide pair, respec- 
tively, A^i, N2 are the ascending nodes, ii and 12 
denotes the inclinations of the orbits measurable 
by an observer while im is the mutual inclination 
of the two orbital planes. Ji and J2 are the peri- 
centre points. 
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Fig. 2. — The UV coverage of the Algol with 
CHARA. Each point represents one observational 
point. 
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Fig. 3. — The result of the grid search for mini- 
mum value of ■ Only the result we obtained for 
f2i is shown here. Note that the logarithm of the 
is shown for sake of a better visualization. 
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Fig. 4. — The modelled versus observed visibility 
values and their errors. The solid line means the 
1:1 line. 
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Fig. 5. — The 6-min. averaged data points with 
our best fit. We also plotted the theoretical posi- 
tions of the points in the case of a revolution in 
a Keplerian orbit in the distance of the secondary 
component. 
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Fig. 7. — The variation of the angular orbital ele- 
ments of Algol system between AD 1600 and 2100. 
The panels from top to bottom: (h) node of the 
binary measured in the invariable plane measured 
from the intersection of the invariable plane and 
the sky; (fi) node of the close binary; (fi') same 
for the third component; (im) mutual inclination; 
(ii) inclination of the binary with respect to the 
invariable plane; (i) inclination of the close binary; 
{i') same for the tertiary. 



Fig. 6. — The astrometric fit with fix a for the data 
given in Lestrade et. al. (1993) (dashed line), as 
well as for the recalculated data, after correction 
of the binary's position along the wide orbit. We 
plotted relative distances with respect to the first 
data point. 
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Fig. 8. — The variation of the inchnation of Algol 
AB with time. Solid line represents the inclini- 
taion of the close pair observable from the Earth. 
Between the dashed lines the system shows par- 
tial eclipses, moreover between the dashed-dotted 
lines the system shows total eclipses for a short 
time. 
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